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ABSTRACT 

We present the analysis of the broadened, flourescent iron Ka line in simultaneous XMM-Newton and RXTE data from the black 
hole Cygnus X- 1 . The XMM-Newton data were taken in a modified version of the Timing Mode of the EPiC-pn camera. In this mode 
the lower energy threshold of the instrument is increased to 2.8 keV to avoid telemetry drop outs due to the brightness of the source, 
while at the same time preserving the signal to noise ratio in the Fe Ka band. We find that the best-fit spectrum consists of the sum 
of an exponentially cut-off power-law and relati vistical 1 y smeared, ionized reflection . The shape of the broadened Fe Ka feature is 
due to strong Compton broadening combined with relativistic broadening. Assuming a standard, thin accretion disk, the black hole is 
close to maximally rotating. 
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1. Introduction 

Relativistically broadened lines from accreting black holes have 
an important diagnostic potential as they allow us to measure 
directly properties of these black holes such as their spin, the in- 
clination of the accretion disk, and the emissivity profile of the 
accretion disk (e.g„ Reynolds & Nowak 2003). While by now a 
large number of high signal to noise ratio (S /TV) measurements 
of the line shape and parameters exist for many Galactic black 
holes (Miller 2007) and for Active Galactic Nuclei (e.g.,Miniutti 
et al. 2007; Fabian et al. 2009), high S /N measurements of the 
line profile and black hole parameters of the canonical black 
hole, Cygnus X-1, are still lacking. This lack is mainly due to 
a combination of source brightness, which causes the measure- 
ment to be dominated by systematic effects in the detector cal- 
ibration (see, e.g.. Miller et al. 2010), as well as a complex X- 
ray spectrum, which necessitates broad band (3-100 keV) mea- 
surements from multiple satellites. See, however. Miller et al. 
(2002) for recent work on the Fe line region with Chandra's grat- 
ings, and Nowak et al. (201 1, and references therein) on joint 
Suzaku/RXTE measurements of the broad band spectrum. 

In this Letter we present an analysis of the broad band spec- 
trum of Cygnus X-1 using data from high S/N observations 
of the Fe line band from XMM-Newton and simultaneous, 3- 
120 keV RXTE data. In Sect. 2 we present the Modified Timing 
Mode, a special flavor of the XMM-Newton Timing Mode, ded- 
icated to bright sources with fluxes up to a few lOOmCrab. In 


Sect. 3 we discuss the data reduction in greater detail. Section 4 
presents the analysis of the relativistically broadened line and 
the reflection continuum. 

2. The EPIC-pn Modified Timing Mode 

To determine the best Fe line parameters possible, very high S /N 
observations are required in a time that is as short as possible to 
avoid complications due to source spectral variability. The best 
instrument currently fulfilling these requirements is the EPIC- 
pn camera onboard XMM-Newton (StrUder et al. 2001). This 
camera has two modes which provide the fast readout capabil- 
ity required to avoid pile up in bright source data (Ness et al. 
2010). Here, the term “pile up” summarizes both "energy pile 
up”, where more than one X-ray photon impacts the same CCD 
pixel during one CCD readout cycle, and "pattern pile up" or 
“grade migration”, where multiple X-ray photons hit adjacent 
CCD pixels during one readout cycle and the resulting charge 
distribution in the CCD is wrongly interpreted as originating 
from a single, higher energy, X-ray. For satellites with moderate 
spatial resolution such as XMM-Newton , the latter effect is more 
pronounced. In EPIC-pn 's Burst Mode, fast exposures and shifts 
are performed, interrupted by longer periods for data read out. 
The livetime of this mode is only 3%. In the Timing Mode, con- 
tinuous read out is available with a 99.5% livetime. This mode 
is well suited for observations with a low pile up fraction for 
fluxes up to ~150mCrab, although the EPIC-pn telemetry allo- 


1 




Duro et al.: Broad Fe K a line of Cygnus X-l 


22:00 23:00 0:00 1:00 2:00 3:00 


1400 

1200 

1000 x 
* 

800 « 

s 

soo jr 

400 
200 
0 

0 5000 10000 15000 

Time |s] 

Fig. 1. XMM-Newton and RXTE lightcurve of Cygnus X-l on 2004 
November 20/21 . Here, and throughout the reminder of the letter, EPIC- 
pn data are shown in blue, PCA (PCU2) data in red, and HEXTE in 
orange, with resolutions of 1 00 s, 96 s, and 40 s, respectively. 



cation limits this mode to observations of sources of fluxes be- 
low ~100mCrab. With ~300mCrab, Cyg X-l can therefore not 
be observed with the standard Timing Mode. 

In order to address the telemetry limitation we have proposed 
a modification of the Timing Mode (Kendziorra et al. 2004). 
In this Modified Timing Mode, the lower energy threshold for 
telemetered events is increased from 0.15 keV to 2.8 keV, such 
that all events in the Fe line band can be transferred without 
telemetry drop outs, and the telemetry allocation to the EPIC-pn 
is maximized by switching off the EPIC-MOS cameras, which 
are piled up anyway for sources above ~35 mCrab. 

EPIC data from events in which the charge cloud produced 
by an X-ray is distributed over more than one CCD pixel, 
so-called “split events”, are only recombined on ground. In 
the Modified Timing Mode, split partners with energies below' 
2.8 keV are not telemetered, resulting in a slight degradation of 
the energy resolution and necessitating the generation of a ded- 
icated response matrix. As discussed by Wilms et al. (2006a) 
and Fritz (2008), this effect can be fully taken into account by 
measuring the energy dependent probability distribution of split 
events from archival Timing Mode observations and modifying 
the official Timing Mode response matrix accordingly. A full de- 
scription of the calibration will be given in a forthcoming pa- 
per (Duro et al. 2011). The XMM-Newton data analysis and 
response matrix generation was based on the Science Analysis 
Software (SAS), version 10.0.0, and the newest calibration files. 

It has recently been found that some Timing Mode observa- 
tions are affected by “X-ray loading” (Done & Diaz Trigo 2009), 
i.e., the contamination of EPIC-pn’s offset map by source X- 
rays. Since the offset map is subtracted from the EPIC data be- 
fore they are telemetered to ground, loading distorts the X-ray 
spectrum and effectively introduces a CCD-column dependent 
lower energy threshold. In the Modified Timing Mode, the off- 
set map is determined with the EPIC-pn filter in the blocking 
position. As no source X-rays enter the offset map. the mode is 
unaffected by X-ray loading. 


3. Data Reduction 

3.1. Data Reduction 

In 2004 fall XMM-Newton performed four Modified Timing 
Mode observations of Cyg X- 1 . Here, we concentrate on the si- 
multaneous XMM-Newton and RXTE observations performed 



Fig. 2. a Fit of a simple power law to the 4.0-9.5 keV data, ignoring the 
5.0-8.0 keV band (light gray area), b A broad excess due to the broad 
line, is visible. Note the energy shift between the PCA and the EPIC-pn. 
c Residuals after applying the gain shift correction. 


2004 November 20 and 21 (XMM-Newton ID 0202760301), 
during which Cyg X-l was in a transitional state between the 
hard state and the thermally dominated soft state (Fritz 2008). 
We concentrate on this 17.4ks long observation since it shows 
the least variability on timescales of minutes to hours (Fig. 1) 
and thus provides the cleanest X-ray spectrum. Results from all 
observations will be presented in a subsequent paper (Duro et al. 
201 1). The increased soft X-ray emission during the intermedi- 
ate state potentially leads to pile-up in the center of the point 
spread function. For the present study, we therefore ignore the 
innermost three columns. In order to avoid any transient effects 
due to the 2.8 keV lower energy threshold, we limit ourselves to 
the 4. 0-9 .5 keV band. 

The simultaneous 6.18 ks of RXTE data (observation IDs 
90104-01-02-00, -01, and -03) were reduced with HEAS0FT 
version 6.9. We use the 4-40 keV data from the Proportional 
Counter Array (PCA; Jahoda et al. 2006) and the 20-112keV 
data from the High Energy X-ray Timing Experiment (HEXTE; 
Rothschild et al. 1998) and ignore data taken within 10 min- 
utes of passages through the South Atlantic Anomaly and dur- 
ing times of high particle background. All spectral analysis was 
done with the Interactive Spectral Interpretation System (Houck 
& Denicola 2000; Noble & Nowak 2008), HEXTE spectra were 
rebinned to S/N = 30, EPIC-pn to S/N = 100, resulting in a 
resolution of 40 eV at 6.4 keV and oversampling the detector re- 
sponse by a factor ~4. 

3.2. Simultaneous XMM-Newton and RXTE Fits 

The 3-120keV spectrum of Cyg X-l can be well described as 
the sum of a soft excess with a temperature of <500 eV, an ex- 
ponentially cutoff power law with a cutoff at a few 1 00 keV, and 
a Compton reflection component (Nowak et al. 2011 , and ref- 
erences therein). Alternatively, a variety of other, more physics 
based, components have been used to describe the continuum 
(Nowak et al. 2011; Maizac et al. 2008; Wilms et al. 2006b; 
Markoff et al, 2005), but these are virtually indistinguishable 
from this empirical model. In addition to the continuum com- 
ponents, Qiandra-HETGS observations show the presence of a 
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narrow line at 6.4 keV emitted by neutral iron in material sur- 
rounding the black hole at larger distances (Hanke et al. 2009). 

As shown in Fig. 2, a simple continuum model is not suf- 
ficient to describe the data. Strong positive residuals in the re- 
gion between 5 and 7 keV remain , which are due to the broad Fe 
K a line. However, a slight mismatch in the RXTE and XMM- 
Newton residuals is also apparent (Figure 2b), with the EPIC 
line appearing to be shifted to a slightly higher energy. This shift 
is confirmed by measuring the energy of the 6.4 keV narrow line 
feature, which is found at 6.6 keV, inconsistent with all available 
Chandra measurements for this feature in Cyg X-l . 

We attribute this difference to the treatment of charge trans- 
fer efficiency (CTE) effects in the SAS. When observing bright 
sources with a CCD such as the EPIC-pn, traps in the Silicon 
crystal can be saturated with source electrons, improving the 
CTE. While the latest SAS versions take this effect into account, 
their calibration breaks down for very high count rates. An over- 
correction of CTE effects then leads to inferred photon energies 
that are too high. We take this effect into account by correcting 
the assigned photon energies with a model of the form £ rea i = 
Eobs/n + A E. The model parameters s and A E are determined 
using measured energies of narrow spectral components or from 
continuum fitting. As will be discussed in greater detail by Duro 
et al. (201 1), Chandra’s HETGS observation (obsid 3814) was 
taken at the same flux level as our XMM-Newton/RXTE obser- 
vations. In addition to the 6.4 keV line, Chandra reveals two Fe 
Ka absorption lines due to Fe xxv and Fe xxvi at 6.646 keV and 
6.955 keV. Fits show that full consistency between the Modified 
Timing Mode spectrum and HETGS is obtained for s = 1.02 and 
A E = 0 keV. This overcorrection of the EPIC energies by 2% is 
in line with typical CTE corrections for bright sources. An initial 
fit to the simultaneous 4-1 OkeV EPIC and PCA data agrees with 
this result (Fig. 2c), showing that the CTE correction can also be 
obtained from continuum fitting. The remaining differences be- 
tween the EPIC and PCA data are due to narrow features which 
are not resolved by the PCA. 

4. The Relativistic Broadened Fe K« Line of Cyg X-1 

We now turn to modeling the broadened Fe line. To describe 
its profile we use the relline-model of Dauser et al. (2010), 
a relativistic line model for thin disks around black holes with 
spins -0.998 < a < +0.998, where a < 0 indicates that the 
angular momenta of the black hole and the disk are antiparallel. 
We assume a radius-dependent line emissivity per disk unit area 
that scales cc r~ e , where e = 3 for a thin accretion disk. The disk 
extends from the (a dependent) innermost stable circular orbit to 
400GM/C 2 . 

As expected, modeling the 4-1 OkeV spectrum with a sim- 
ple relativistic line and a power-law' continuum results in a bad 
fit fyjC ~ 1 .5) and unphysical line parameters as this approach 
does not take reflection self-consistently into account. We there- 
fore extend the data to the full 4— 120keV band provided by 
EPIC-pn, PCA, and HEXTE, and model reflection of an ex- 
ponentially cutoff power law self-consistently using the reflec- 
tion model reflionx (Ross & Fabian 2005). This model in- 
cludes transitions from the most important ions, including Fevt 
to Fexxvi. In order to account for the relativistic smearing we 
convolve the reflected spectrum with a relativistic convolution 
model (relconv: Dauser et al. 2010). We assume that the an- 
gular distribution of the fluorescent photons is isotropic and due 
to photons emerging from a hot corona (Svoboda et al. 2009). 
In addition, we include narrow Ka absorption lines from Fexxv 
and Fexxvi and a soft excess in the model. 
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Fig.3. a Unfolded XMM-Newton and RXTE data and best fit model 
components. Green line: reflection continuum in the frame of the disk 
(note the strong Compton broadening; the “jitter’ above 50 keV is due 
to the numerical resolution of the reflection model). Purple line: rela- 
tivistically smeared reflection component, b Measured count rate spec- 
tra. c Residuals of the best fit setting the relativistic convolution to zero, 
d Best fit residuals including relativistic convolution. 

We first model the continuum using the canonical value of 
e = 3 for a thin accretion disk. As shown in Fig. 3 and Table 1 , 
this model describes the 4— 120 keV energy spectrum well, with 
a reduced x^ cii = 1.10. The line profile indicates a high angu- 
lar momentum of the black hole ( a = 0.88*q®, Fig. 4b), with 
further line broadening being due to strong Compton broaden- 
ing in the strongly ionized reflector (ionization parameter £ = 
1 400C[Jj erg cm s _1 ) , which smears out all other discrete features 
such as edges in the reflection spectrum (Fig. 3, green line). The 
inferred inclination is 32° ± 2°, in good agreement with most es- 
timates for the inclination of the system (e.g., 32° < i < 40°, 
Ninkovetal. 1987,1 < 55°, Sowers et al. 1998, or 36° < i < 67° 
Davis & Hartmann 1983. The Fe abundance is slightly higher 
than solar. 

In order to check the uniqueness of this thin disk solution, 
a second fit letting e free was performed. Figure 4 shows sig- 
nificance contours for this case. Most best-fit parameters are 
consistent between both fits, specifically, the inclination remains 
around 35° and the reflector remains strongly ionized. The emis- 
sivity, however, increases to e = 10 (the hard limit) while a 
decreases to 0. Broad lines can therefore be generated either 
by strongly concentrating all available emissivity to the inner- 
most regions of a disk around a Schwarzschild black hole, or 
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Table I. Best fit results to gainshifted F.PIC-pn, PCA, and HEXTE 
data assuming an exponentially cutoff power law, a narrow emission 
line at 6.4 keV, two resonant absorption lines, disk black body emis- 
sion, and relativisticaly convolved reflection. Fluxes and normaliza- 
tions, A,., are determined relative to the RXTE-PCA. Cehc and Chexte 
are flux normalization constants for the EPIC and HEXTE, respectively. 
Uncertainties are at the 90% level for one interesting parameter. 
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e frozen 
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J K 534 -O.OD 
1 ‘ WJ -0.0|7 
1 17+0.09 
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290/f 

290xf 

Abb 

(0.44:™) xl(F 
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E(.A kcY [cgs] 
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(0.69!™) x 10" 3 

Ep e XXV Kir [KcVJ 
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T l 
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(1.4 + 0.3) x 10“ 

j 7+0,5 

Mr— 
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-0.1 ±0.4 
36 + j 


32 ±2 


0-830:^ 
0.785 ± 0.004 
1 024O +00019 


Chexte 

cepic 

■^Kainshift 


X^/dof 

Knell 


0.830 ± 0.005 
0.785 ± 0.004 
l-0230^_ 


254/237 

1.08 


261/238 
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EPIC-pn camera to observe sources up to several lOOmCrab, 
yielding high S/N measurements of the Fe line band. Coupled 
with broad band data from other satellites such as RXTE, 
INTEGRAL, or Suzaku, the Modified Timing Mode allows the 
high precision measurement of the physical parameters of black 
hole candidates, including their inclination, spin, and emissivity 
profile. Applying this analysis to Cyg X-l shows that its black 
hole is likely a Kerr black hole if one assumes the accretion disk 
to be flat and illuminated from a Compton corona. This result is 
consistent with spin measurements from the accretion disk con- 
tinuum (Gou et al. 201 1). Being in a High Mass X-ray Binary, 
the black hole in Cyg X-l is comparatively young. Our result 
indicates that it must have been born at close to maximum spin 
as expected in typical scenarios for black hole formation in su- 
pernovae (MacFadyen & Woosley 1999). In a next step, all four 
Modified Timing Mode observations will be considered (Duro 
et al. 201 1). 
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Fig. 4. a x 2 significance contours for the emissivity index, e, and the 
spin, a. Contours are based on Ax 2 = 2.30, 4.61, and 9.21, i.e., 68%, 
90%, and 99% confidence for two parameters of interest, with respect to 
the best fit case with e let free. Two minima are apparent, one indicating 
a low spin black hole at an unphysically large e, and one indicating an 
almost maximally spinning black hole consistent with emission from a 
thin accretion disk (e = 3). bx 2 as a function of a for e = 3. 


by emission from a standard disk around a maximally rotating 
black hole. Both solutions are comparable in statistical quality 1 . 
Given that most physical scenarios for very steep emissivity pro- 
files such as strongly torqued accretion disks (Agol & Krolik 
2000) or strong illumination of the innermost disk due to strong 
light bending from a high latitude source (“lamppost models”; 
Martocchia & Matt 1996) also require high a , the high spin so- 
lution is preferred on physical grounds. Cygnus X-l is therefore 
plausibly close to maximally rotating. 

To conclude, this first analysis of a data set taken in the 
EPIC-pn Modified Timing Mode confirms the capability of the 

1 For Ax 2 = 7 and 235, respectively 236 degrees of freedom, the F- 
test indicates no improvement by introducing a free disk emissivity at 
the 99% level. 
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